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Abstract. We study the transition from a prestellar core to a Class protostar, using SPH to simulate the dynamical evolution, 
and a Monte Carlo radiative transfer code to generate the SED and isophotal maps. For a prestellar core illuminated by the 
standard interstellar radiation field, the luminosity is low and the SED peaks at ~ 190 yum. Once a protostar has formed, the 
luminosity rises (due to a growing contribution from accretion onto the protostar) and the peak of the SED shifts to shorter 
wavelengths (80 to 100 pm). However, by the end of the Class phase, the accretion rate is falling, the luminosity has decreased, 
and the peak of the SED shifts back towards longer wavelengths (90 to 150/um). In our simulations, the density of material 
around the protostar remains sufficiently high well into the Class phase that the protostar only becomes visible in the NIR 
if it is displaced from the centre dynamically. Raw submm/mm maps of Class protostars tend to be much more centrally 
condensed than those of prestellar cores. However, when convolved with a typical telescope beam, the difference in central 
concentration is less marked, although the Class protostars appear more circular. Our results suggest that, if a core is deemed 
to be prestellar on the basis of having no associated IRAS source, no cm radio emission, and no outflow, but it has a circular 
appearance and an SED which peaks at wavelengths below ~ 170yum, it may well contain a very young Class protostar. 

Key words. Stars: formation - ISM: clouds-structure-dust - Methods: numerical - Radiative transfer - Hydrodynamics 



1. Introduction 

The details of how molecular clouds form, why they collapse, 
and how the collapse proceeds to form stars are not very well 
understood. Over the last two decades, observations of the dif- 
ferent stages of this process have lead to an evolutionary sce- 
nario of star formation: starless core — > prestellar core — > class 
object — > class I object — > class II object — > class III object 
(see Andre et al. 2000). 

Starless cores (e.g. Myers et al. 1983; Myers & Benson 
1983) are dense cores in molecular clouds in which there is 
no evidence that star formation has occurred (i.e. no IRAS de- 
tection, Beichman et al. 1986). Some of these starless cores 
are thought to be close to collapse or already collapsing, and 
they are labelled prestellar cores (Ward-Thompson et al. 1994, 
1999); for the rest, it is not always clear whether they are in hy- 
drostatic equilibrium (e.g. Alves et al. 2001 or transient struc- 
tures within a turbulent cloud (e.g. Ballesteros-Paredes et al. 
2003). Prestellar cores have been observed both in relative iso- 
lation (e.g. L1544, L63; Ward-Thompson et al. 1999) and in 
protoclusters (e.g. p Oph, Motte et al. 1998; NGC 2068/2071, 
Motte et al. 2001). They have typical sizes (1 -20)xl0 3 AU and 
masses from 0.05 to 35 Mq. From statistical arguments (e.g. 
Andre et al. 2000) it is inferred that prestellar cores live only a 
few million years. 
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Class objects represent the stage when a protostar has 
just been formed in the centre of the core, but the protostar is 
still less massive than its surrounding envelope (Andre et al. 
1993). The protostar is deeply embedded in the core and can- 
not be observed directly, but its presence can often be inferred 
from bipolar molecular outflows or compact centimetre radio 
emission. The Class phase last for a few 10 4 years and is 
characterised by a large accretion rate (^ 10~ 5 MQ/yr). Most 
of the mass of the protostar is delivered during this phase (e.g. 
Whitworth & Ward-Thompson 2001). 

Once the protostar becomes more massive than its sur- 
rounding envelope, the core enters the Class I phase. Accretion 
onto the central protostar, or onto the disc around the protostar, 
continues, but at a lower rate (^ 10~ 6 M0/yr). This phase lasts 
for a few 10 5 years, and delivers most of the rest of the mass of 
the protostar. 

Finally, Class II and Class III objects correspond to even 
later stages of star formation, when most of the envelope has 
disappeared, and the star is visible in the optical. Class I ob- 
jects are classical T Tauri stars (CTTSs), with ongoing accre- 
tion (but at much lower levels than before, ^ 10~ 8 MQ/yr) and 
well defined discs. Class III objects are weak-line T Tauri stars 
(WTTs) with little sign of accretion and no inner disc. 

In their quest to identify the youngest protostars, Andre 
et al. (1993, 2000) set 3 criteria for classifying an object as 
Class 0: (i) the presence of a central luminosity source, as indi- 
cated by the detection of a compact centimetre radio source, or 
a bipolar molecular outflow, or internal heating; (ii) centrally 
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peaked but extended submillimetre continuum emission, indi- 
cating the presence of an envelope around the central source; 
and (iii) a high ratio (> 0.005) of submillimetre to bolometric 
luminosity, where the submillimetre is defined as A ^ 350 /mi. 
Criteria (ii) and (iii) are used to distinguish between Class 
and Class I objects, whereas criterion (i) is used to distinguish 
between Class objects and prestellar cores. However, crite- 
rion (i) may be inadequate to distinguish between the youngest 
Class objects and the oldest prestellar cores: (a) the cm emis- 
sion may be be too weak to be observed, especially in the ear- 
liest stages of protostar formation; (b) the region may be too 
complex for the molecular outflows to be detectable; and (c) 
internal heating may be difficult to establish, if the protostar is 
deeply embedded in the core. 

The origin of radio emission from Class objects is un- 
certain (cf. Gibb 1999). Theoretical models show that once a 
protostar forms at the centre of a collapsing core, the infalling 
gas accelerates to supersonic velocities and an accretion shock 
develops on the surface of the protostar and/or on the surface 
of the disc that surrounds the protostar (Winkler & Newman 
1980, Cassen & Moosman 1981). The heating provided by 
the shock ionises the surrounding gas, which then emits free- 
free radio radiation (e.g. Bertout 1983; Neufeld & Hollenbach 
1994, 1996). The radio luminosity depends on the mass of the 
protostar, the accretion rate and the observer's viewing angle. 
For a low-mass protostar, $ Mq, the predicted flux is below de- 
tection limits (Neufeld & Hollenbach 1996), unless the accre- 
tion rate is high, ^ 10~ 5 Mo/yr. Another possibility is that the 
radio emission is produced by an ionised disc wind (e.g. Martin 
1996), or from a partially ionised jet that emanates from the 
protostar and propagates into the collapsing envelope produc- 
ing shocks (Curiel et al. 1987). The detection of radio emission 
from many Class protostars, using the VLA (e.g. Bontemps 
et al. 1995), points toward the latter two explanations. 

Recently Young et al. (2004), using the Spitzer space tele- 
scope, detected NIR radiation from L1014, a dense core that 
was previously classified as prestellar. The detection of NIR 
radiation strongly suggests that this is a young Class object, 
and supports the view that even younger Class objects may 
not be detectable in the NIR, with IRAS, or even with Spitzer. 

The goal of this paper is to study the transition from prestel- 
lar cores to Class objects, and to seek new criteria for distin- 
guishing between these two stages, in particular criteria which 
can be used before the standard signatures of protostar birth, 
such as compact radio emission and/or bipolar molecular out- 
flows, become detectable. We use SPH to simulate to the dy- 
namics of a collapsing molecular core, and a Monte Carlo code 
to treat the transfer of radiation within the core. The radiative 
transfer is treated fully in 3D, which is important for the cor- 
rect interpretation of the observations (cf. Boss & Yorke 1990; 
Whitney et al. 2003a, 2003b; Steinacker et al. 2004). We use a 
newly developed method for performing Monte Carlo radiative 
transfer simulations on SPH density fields, which constructs ra- 
diative transfer cells using the SPH tree structure (Stamatellos 
2003; Stamatellos & Whitworth 2005, hereafter Paper I). 

In Section 2 we describe the SPH simulation of the col- 
lapse of a turbulent molecular core. In Section 3 we describe 
the radiative transfer method, focusing on how we construct 



radiative transfer cells. We also discuss the radiation sources 
and the properties of the dust used in our model. In Section 
4 we present our results for the dust temperature fields, SEDs 
and isophotal maps of prestellar cores and young Class pro- 
tostars. In Section 5 we discuss how SEDs and isophotal maps 
at submm and mm wavelengths might be used to distinguish 
between late-phase prestellar cores and early-phase Class ob- 
jects. In Section 6 we summarise our results. 

2. The SPH simulation 

We use the Smoothed Particle Hydrodynamics code dragon 
(Goodwin et al. 2004a) to simulate the collapse of a turbu- 
lent molecular core and the resulting star formation. Here, we 
briefly describe the main elements of the model. For a more 
detailed discussion we refer to Goodwin et al. (2004b). 

The initial conditions in the core, before the start of col- 
lapse, are dictated by observations of prestellar cores (e.g. 
Ward-Thompson et al. 1999, Alves et al. 2001). Prestellar cores 
appear to have approximately uniform density in their central 
regions, and the density then falls off in the envelope. If the 
density in the envelope is fitted with a power law, n(r) <x r~ n , 
then j] ~ 2 to 4 . Here r\ ~ 2 is characteristic of more extended 
prestellar cores in dispersed star formation regions (e.g. L1544, 
L63 and L43), whereas r\ ~ 4 is characteristic of more compact 
cores in protoclusters (e.g. p Oph and NGC2068/2071). These 
features can be represented by a Plummer-like density profile 
(Plummer 1915), 



The density profile is almost flat for r < r , and falls off as r 
in the outer envelope (r > r ). We set p a — 3 x 10~ 18 gcirT 3 , 
r = 5, 000 AU, and r B = 50, 000 AU. The core mass is then 5.4 
Mq. The core is initially isothermal at 10 K, hence the ther- 
mal virial ratio is a mESM = ^therm/I^gravI = 0.45. Observations 
of star forming cores (e.g. Jijina et al. 1999) show that their 
line widths have a non-thermal contribution, which may be due 
to turbulence. We impose a low level of turbulence (ff TURB = 
E rmB /\E 0RAV \ = 0.05) with a power spectrum P(k) oc k~ 4 so as 
to mimic the observed scaling relation between linear size and 
internal velocity dispersion (Larson 1981). 

Collapse of the turbulent core leads to the formation of a 
single star surrounded by a disc, and the material in the disc 
then slowly accretes onto the star. A sink particle, represent- 
ing the star and the inner part of the disc, is created where the 
density first exceeds p CRIT = 10~" g cirT 3 , and incorporating all 
the matter within R mn = 5 AU of this point (Bate et al. 1995). 
(The values of p CRIT and R mlT are dictated by computational con- 
straints. If p CRIT is increased, R lmT must be reduced, and more 
computational time is needed.) 

We use a newly developed type of sink called a smartie, 
which is a rotating oblate spheroid with radius R s and half- 
thickness Z s (Fig.Q. The star is a point mass M t at the centre 
of the smartie, and the spheroid represents the inner part of the 
accretion disc surrounding the star. R s and Z s are calculated by 
assuming that the inner part of the accretion disc (the part in- 
side the spheroid) has uniform density, temperature and angular 
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Fig. 1. We use an oblate spheroidal sink particle, called a smar- 
tie, which contains a central star and a disc, and launches a 
bipolar outflow. 

speed, and that its internal energy and spin angular momentum 
balance its self-gravity and the gravity of the central star. 

The SPH particles which enter the smartie during a 
timestep, (t, t + Af), are assimilated by it at the end of the 
timestep, t + Af, and their mass is initially added to the mass 
of the smartie, M Q , i.e. to the mass of the inner disc, 



M' = M c 



rtii) 



(2) 



At the same time, the centre of mass, R s , and velocity, V s , of 
the smartie are adjusted to conserve linear momentum, 



R' 



V' 



M S R S + £ ; {ro f r,-) 
M' 

S 

M' 



(3) 
(4) 



and the intrinsic angular momentum (spin) of the smartie is 
adjusted to conserve angular momentum, 



h; = h s + Ms (r s -r;)®(v s -v;) 
+ ^{,«,(r,- R ;)®(v,-v;)}. 



(5) 



In Eqns, |3 through |5] primed variables (M' s , R' s , V', H' s ) are 
used to denote sink parameters adjusted for the addition of as- 
similated SPH particles. 

The smartie mass evolves according to 



dt 



dM 
dt 



dM^ 
dt 



dM 
dt 



(6) 



where the first term on the right of Eqn. ^represents SPH par- 
ticles assimilated by the smartie, as discussed in the preceding 
paragraph. 

The second term on the right of Eqn.[6]represents accretion 
onto the central star, and is given by 



dM t 
dt 



0.9M, 



[G(M+ +M s )R t 



lV2 



1 + 



GM 2 



(M t + M S )R a 2 



(7) 
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where a ss is the standard Shakura & Sunyaev (1973) turbu- 
lent viscosity parameter and a s is the mean sound speed in the 
smartie. The leading term represents turbulent viscosity, and 



the term in braces adjusts this for angular momentum transport 
by gravitational torques in a Toomre-unstable disc (cf. Lin & 
Pringle 1990). 

The third term on the right of Eqn. ^represents mass loss 
from the smartie, which can occur in two ways, (i) SPH parti- 
cles are created at rate 0.1M s /f VISC and launched along the ro- 
tation axis with speed 100kms~' to form a bipolar outflow. 
These jets push the surrounding material aside and create hour- 
glass cavities about the rotation axis, (ii) In the present model 
the outflow does not remove any angular momentum. Therefore 
the redistribution of angular momentum which drives accretion 
onto the central star causes the residual inner disc (which has to 
carry all the angular momentum assimilated by the smartie) to 
expand. When this happens, the smartie is allowed to expand 
until R s = 16 AU. Thereafter, R s is held constant by excret- 
ing SPH particles - and with them angular momentum - just 
outside the equator of the smartie. 

The luminosity of the star is given by 



GM t dM t 
~R ~dT 



(9) 



where M t and R t are the mass and radius of the star, and we 
set/?, = 3R Q (Sta'hler 1988). 

The mean smartie temperature, and hence also the mean 
smartie sound speed, are determined by the balance between 
radiative cooling, heating by radiation from the star and from 
the background radiation field, and viscous dissipation in the 
inspiralling gas of the smartie. A full description of the imple- 
mentation of smarties is given in Boyd (2003), and the con- 
sequences of feedback from bipolar outflows are explored in 
detail in Boyd et al. (in preparation). 

3. The radiative transfer simulation 

For the radiative transfer simulations, we use a version of 
phaethon, a Monte Carlo radiative transfer code which we have 
developed (Stamatellos & Whitworth 2003) and optimised for 
radiative transfer simulations on SPH density fields (Paper I). 
The main developments of the optimised version are (i) that 
it uses the tree structure inherent within the SPH code to con- 
struct a grid of cubic radiative transfer (RT) cells, spanning the 
entire computational domain, the global grid; and (ii) that in 
addition it constructs a local grid of concentric spherical RT 
cells around each star, a star grid, to capture the steep tem- 
perature gradients that are expected in the vicinity of a star. 
The code reemits luminosity packets are soon as they are ab- 
sorbed (Lucy 1999) and uses the frequency distribution adjust- 
ment technique developed by Bjorkman & Wood (2000). 

3.1. Construction of the global radiative transfer grid 

To construct RT cells for a given time-frame from an SPH sim- 
ulation, we take advantage of the fact that SPH uses an oc- 
tal tree structure (to group particles together when calculat- 
ing gravity forces and also to find neighbours; for details see 
Paper I). The SPH tree is a recursive hierarchical division of 
the computational domain into cubic cells within cubic cells 
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regular radiative transfer 
cells 




star grid 



Fig. 2. The star grid consists of concentric spherical surfaces 
centred on the star, with equal spacing in log(r), and conical 
surfaces symmetric about the angular momentum of the smar- 
tie, with equally spaced opening angle 6. This grid takes ac- 
count of the dust-free cavity in the immediate vicinity of the 
star, and resolves the disc inside the smartie. 



(Barnes & Hut 1986). When the SPH tree is being built, we 
record information about the size of each cell and the number 
of SPH particles it contains. The RT cells of the global grid are 
then the largest SPH cells which contain < N MAX particles. This 
means that the RT cells automatically have comparable mass. 
We choose A^ MAX ~ 100 to 150, i.e. somewhat larger than the 
mean number of SPH neighbours, N NEm ~ 50. Consequently 
the size of the RT cells is on the order of the smoothing length, 
h, and so the temperature resolution is similar to the density 
resolution of the SPH simulation. 

In order to increase the resolution of the radiative transfer 
simulation, we use particle splitting (Kitsionas & Whitworth 
2002). This method replaces each SPH particle (parent particle) 
with a small group of 13 particles (children particles), each one 
having 1/13 of the mass of the parent particle. The children are 
placed on an hexagonal close-packed array, with one of them 
in the centre of the array and the other 12 equidistant from the 
first one. By doing this the smoothing length is decreased by a 
factor of 13~ 1/3 = 0.425, and the average size of the radiative 
transfer cells is reduced by the same factor. 

3.2. Construction of the local radiative transfer grid 
around a star 

In Paper I we discussed the need for a high-resolution grid of 
radiative transfer cells - a star grid - in order to capture the 
steep temperature gradients that are expected in the vicinity of 
a star. The star grid discussed there was developed for treating 
the unresolved region inside a spherically symmetric sink, and 
so a one-dimensional star grid consisting of concentric spheri- 
cal cells centred on the star was sufficient. 

In the case we examine here the sinks are smarties, i.e. 
oblate spheroids, and the unresolved region inside a smartie 
represents both the protostar and the inner protostellar disc. 



Thus a spherically symmetric one-dimensional star grid is in- 
appropriate. Instead we adopt a similar approach to Kurosawa 
et al. (2004) and represent the unresolved region near the star 
with a flared disc (Fig.|5Jl. 

However, the properties of the inner disc are not imposed 
externally, as in Kurosawa et al. (2004), but are determined by 
the properties of the smartie, and therefore attempt to capture 
the physics of a protostellar disc accreting onto a newly-formed 
protostar (as outlined in Section Specifically, we adopt an 
inner disc with inner radius R^ T , outer radius /?„ , and constant 

DEST 1 SG 1 

ratio of thickness, H, to radius, r, 



H(r) 



a r , where a 



5l 



The dust destruction radius is given by 



D _ * 



(4+«/2 



(10) 



(11) 



where /3 is the dust emissivity index (see Paper I). There is no 
dust inside R VEST . 

We put Rsg = 16 AU, since this is the maximum size of a 
smartie, and it is comparable with the size of the local global 
RT cells (the cubic cells constructed from the SPH tree). 

For z < 4ar and ^ DEST < r < R sa , the density of the inner 
disc is given by 



Piz, r) = p D 



exp 



2 Karl 



(12) 



(cf. Wood et al. 2002a, 2002b). p DEST is fixed by equating the 
mass of this inner disc to the mass of the disc inside the smartie. 
For z > 4-ar and R DEST < r < R SG , the density is set equal to the 
average density of the local global RT cells. 

The inner disc is then divided into RT cells by spherical and 
conical surfaces (see Fig. [2}. The spherical surfaces are evenly 
spaced in logarithmic radius and there are typically 20 of them. 
The conical surfaces are evenly spaced in polar angle and there 
are typically 30 of them. 

3.3. Radiation sources 

The core is illuminated externally by the interstellar radiation 
field and internally by the newly formed protostar (once it has 
formed). 

For the external radiation field we adopt a revised version 
of the Black (1994) interstellar radiation field (BISRF). The 
BISRF consists of radiation from giant stars and dwarfs, ther- 
mal emission from dust grains, cosmic background radiation, 
and mid-infrared emission from transiently heated small PAH 
grains (Andre et al. 2003), as illustrated on Fig. [3] Typically the 
BISRF is represented by 10 9 - 10 10 L-packets (see Paper I). 

The protostar luminosity, L^, and effective temperature, 
T t are given by Eqn. |9j a blackbody spectrum at T t is as- 
sumed. For the systems examined here, the accretion con- 
tribution to the luminosity dominates, because the mass of 
the protostar is small (< 0.6 Mq) and the accretion rate is 
high (> 10~ 7 Mo/yr). Typically the protostar is represented by 
10 6 to 10 8 L-packets. 
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Table 1. Model parameters 



id 


time (yr) 


M,(M Q ) 


As (AU) 


Z s (AU) 


M S (M ) 


M,(M /yr) 


^tot(L©) 


T* (K) 


description 


to 



















initial conditions 


tl 


5.3 x 10 4 
















before sink formation 


t2 


5.3 x 10 4 


0.01 
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3.4 


0.04 


5 x 10~ 5 


5.7 


5150 


a protostar has formed 


t3 


6.0 x 10 4 


0.20 


4.0 


0.4 


0.09 


1 x 10~ 5 


27.2 


7650 




t4 


6.8 x 10 4 


0.49 


4.0 


0.4 


0.05 


2 x 10- 6 


12.3 


6250 




ts 


6.9 x 10 4 


0.53 


4.0 


0.4 


0.01 


4 x 10~ 7 
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4200 


protostar is off centre 
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Fig. 3. Dotted line: The external illumination field (BISRF + 
emission from PAHs). Solid line: The radiation scattered by 
the molecular cloud. This is similar for all the time-frames and 
viewing angles, for the models we examine here. 



3.4. Dust opacity 

Like other studies of prestellar cores (e.g. Evans et al. 2001, 
Young et al. 2004) we use the Ossenkopf & Henning (1994) 
opacity for a standard MRN interstellar grain mixture (53% sil- 
icate and 47% graphite) which has coagulated and accreted thin 
ice mantles over a period of 10 5 yr at a density of 10 6 cm -3 . 
However, we emphasize that the opacity of the dust in cores is 
very uncertain (e.g. Bianchi et al. 2003). 



4. The evolution of the dust temperature field, SED 
and isophotal maps in a star-forming core 

We perform radiative transfer simulations on 6 time-frames 
during the collapse of a star-forming core. We focus our atten- 
tion just before and just after the formation of the first protostar 
in the core. The characteristics of each time-frame are listed in 
Tabled The simulations are 3-dimensional and provide dust 
temperature fields, SEDs and isophotal maps. 



4.1. Dust temperature fields in prestellar cores and 
young protostars 

In Fig. |4] we present the density fields and dust- temperature 
fields for the 6 time-frames in Tabled We plot the density and 
the temperature on 2 cuts through the computational domain, 
one parallel to the x = plane and one parallel to the z = 
plane. Each of these planes passes through the protostar, or, if 
there is no protostar (as in the first 2 time-frames), through the 
densest part of the core. The plots show only the central region 
of the core (5000 AU x 5000 AU). Additionally, on Fig. [5] we 
plot the density and the temperature of each RT cell versus dis- 
tance r from the centre of coordinates, on a logarithmic scale, 
in order to depict better the regions very close to the protostar. 

Our results for the temperature are broadly similar to those 
of previous ID and 2D studies of prestellar cores and Class 
objects. Before the collapse, the precollapse prestellar core is 
quite cold (10 to20K) and it becomes even colder (5 to20K) 
as the collapse proceeds and the central regions become very 
dense and opaque. As soon as a protostar forms, the region 
around it becomes very hot (up to the dust destruction tem- 
perature), but the temperature still drops below ~ 20 K beyond 
a few 1000 AU from the protostar, because of the high optical 
depth in the dense accretion flow onto the protostar. The lu- 
minosity of the protostar is dominated by the contribution from 
accretion. Initially the accretion luminosity increases due to the 
increasing mass of the protostar, but then it falls as the accre- 
tion rate declines. 

4.2. SEDs of prestellar cores and young protostars 

The SEDs of the 6 time-frames in Table ^ are presented in 
Fig. |6] These SEDs have been calculated assuming that the 
core is at 140 pc. SEDs are plotted for 6 different viewing an- 
gles, i.e. 3 polar angles (6 = 0°, 45°, 90°) and 2 azimuthal 
angles (<p = 0°, 90°). Hence on each graph there are 6 curves. 
However, in several cases the SED is so weakly dependent on 
viewing angle that the curves overlap. We plot only the ther- 
mal emission from the core, i.e. we neglect both the long- 
wavelength background radiation that just passes through the 
cloud, and the short-wavelength radiation that is scattered by 
the outer layers of the cloud. 

The effective temperature of the core, as inferred from the 
peak of the SED, rises and falls with the accretion luminosity of 
the protostar. For a prestellar core, the SED peaks at ~ 190 pm, 
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Fig. 4. Cross sections of density and dust temperature on a plane parallel to the x — plane (left two columns) and parallel to 
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implying T EFF ~ 13 K. Once the protostar has formed and in- 
puts energy into the system, the peak moves steadily to shorter 
wavelengths, reaching ~ 80 /mi (!T EFF ~ 3 IK) as the accre- 
tion luminosity reaches its maximum, and then moving back to 
longer wavelengths again as the accretion luminosity declines. 
By the final frame it has reached ~ 150 //m (r EFF ~ 17 K). 

The peak of the SED of a prestellar core is independent of 
viewing angle, since the core is optically thin to the radiation it 
emits. In contrast, the peak of the SED of a Class object does 
depend on viewing angle, albeit weakly, because of the pres- 
ence of an optically thick disc around the protostar. However, 
even allowing for variations in the viewing angle, the SED of a 
Class object does not peak at the wavelengths characteristic 
of prestellar cores (~ 190 //m). 

4.3. Is a young Class observable in the NIR? 

Our model predicts that a young protostar embedded in a core 
is not observable in the NIR, unless the protostar is displaced 
from the central high-density region. This contradicts the recent 
results of Whitney et al. (2003b) and Young et al. (2004). We 
attribute the difference between our model and that of Young 
et al. to the density they use for the region around the proto- 
star (~ 10 5 cirT 3 ), which may be too low for a Class ob- 
jects and more appropriate for a prestellar core. Whitney et 
al. use a higher density in the central region (~ 10 10 cirT 3 ), 
but it is still more than one order of magnitude lower than the 
density produced by our SPH simulation (~ 3 x 10 12 cirT 3 ). 
Additionally, they use different dust properties for the differ- 
ent regions in their configuration (cloud, disc mid-plane, disc 
atmosphere, outflow); in particular, the opacity of their dust in 
the disc midplane is smaller than the opacity in our model, at 
wavelengths shorter than 10 /mi. There are also some differ- 
ences in the properties of the discs used in the different models, 
but these are less relevant to the escape of NIR radiation. 

The models reported by Shirley et al. (2002) use central 
densities and dust opacities similar to our models, but they con- 
fine their study to the FIR/mm region of the spectrum, so a 
comparison is not possible. 

Ultimately, the density distribution within 500 AU of the 
centre of a core or protostar is not well constrained by observa- 
tions, and therefore radiative transfer calculations have to rely 
on theoretical models of core collapse. Dust opacities are also 
poorly constrained (e.g. Bianchi et al. 2003). Comparison of 
our model with those of Whitney et al. (2003b) and Young et 
al. (2004) suggests that the presence or absence of NIR emis- 
sion from a confirmed Class object might be used to constrain 
the density profile within a few 100 AU of the protostar and/or 
the opacity of the dust in the same region. 

4.4. Isophotal maps of prestellar cores and young 
protostars 

The shapes of prestellar cores and Class objects depend on 
viewing angle. In Fig.Qwe present isophotal maps at 850 /mi, 
for 3 different time-frames (tl, t3, t5), and from 3 different 
viewing angles. At 850 /mi the core is optically thin, and the 
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Fig. 6. SEDs of the 6 time-frames in Tabled from 3 polar angles (8 = 0°, 45°, 90°) and 2 azimuthal angles {<p = 0°, 90°); in 
several cases the SEDs from different viewing angles overlap, (a) tO - precollapse prestellar core. The SED peaks at ~ 190 pm, 
which implies an effective temperature T EEE ~ 13 K. (b) tl - collapsing prestellar core. The SED appears the same as in (a) 
although the collapse has started and the cloud is denser and colder at its centre, (c) t2 - Class object. The emission of the 
system peaks at A ~ 100 to 135ji/m (depending on the observer's viewing angle), which implies r EFF ~ 18to25K. The total 
luminosity of the system is higher than in (b). (d) t3 - Class object. Similar to (c), but the luminosity is even higher, and the 
SED peaks at A ~ 80 to 100/im, implying T EEE ~ 25 to 31 K. (e) t4 - Class object. Similar to (d), but the luminosity has started 
to decrease due to the falling accretion rate, and the SED peaks at A ~ 90 to 135 pm, implying T EEE ~ 18 to 28 K. (f) t5 - Class 
object. The protostar has moved out of the central region by ~ 100 AU, and the attenuated stellar emission can be seen at short 
wavelengths (1 to 50 pm). The peak of the cloud emission is at A ~ 150 pm, implying T EEE ~ 17 K. 



temperature does not vary much (~ 10 to20K, apart from the 
region very close to the protostar in time-frames t3 and t5), so 
the maps are effectively column density maps. Class objects 
are more centrally condensed than prestellar cores. They also 
show more structure, due to bipolar outflows, which clear low- 
density cavities (see images at — 90° and = 90° on Fig. 0. 

In Fig. [8] we present a selection of isophotal maps af- 
ter convolving them with a Gaussian beam having FWHM = 
2,300AU. Assuming the cloud is 150 pc away, this corre- 
sponds to an angular resolution of 15", which is similar to the 
beam size of current submm and mm telescopes (e.g. SCUBA, 
IRAM). On these maps, Class objects look very similar to 
prestellar cores, apart from the fact that they tend to appear 
more circular and featureless (cf. Fig. 0. 

The reason for this is that the emission from a core that 
contains a protostar is concentrated in the central few hundred 
AU, and so, when it is convolved with a 2, 300 AU beam, it 
produces an image rather like the beam, i.e. round and smooth. 



In contrast, the emission from a prestellar cores has structure 
on scales of several thousand AU, much of which survives con- 
volution with a 2, 300 AU beam. Thus we should expect cores 
that contain protostars to appear rounder than prestellar cores, 
and indeed this is what is observed (Goodwin et al. 2002). 

4.5. Individual time-frames 

Time-frame tO, precollapse prestellar core. The core is heated 
exclusively by the ambient radiation field, and so the outer parts 
of the cloud are warmer (~ 20 K) than the interior. The tempera- 
ture falls to ~ 8 K in the central ~ 5000 AU of the cloud, where 
the density is high, ~ 10 6 ctrT 3 (cf. Evans et al. 2001, Zucconi 
et al. 2001, Stamatellos & Whitworth 2003, Goncalves et al. 
2004, Stamatellos et al. 2004). The cloud emission peaks at 
~ 190 pm which implies 7^ ~ 13 K. The SED is fairly typical 
for a prestellar core (e.g. Ward-Thompson et al. 2002), except 
that T EEE is a little higher than for the majority of prestellar cores 
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(where it is ~ 10 K). This is because the SED plotted here in- 
cludes the outer, warmer layers of the core, and this shifts the 
peak of the SED to shorter wavelengths. 

Time-frame tl, collapsing prestellar core. The collapse of 
the core results in a flattened, disc-like region in the centre of 
the cloud. Heating is still provided only by the ambient radia- 
tion field. The temperature at the centre of cloud is even lower 
(~ 5 K) than in time-frame tO, because the density - and hence 
the optical depth - is even higher (> 10 s cm 3 ). The SED looks 
very similar to time-frame tO, even though the collapse has 
started and the cloud is denser and colder at its centre. Thus the 
SED does not distinguish between a precollapse prestellar core 
and a collapsing prestellar core. 

Time-frame t2, Class object. A protostar has formed at 
the centre of the core, and started heating the core. The proto- 
star has low mass (~ 0.01 Mq), but the accretion rate onto it is 
high (M* ~ 5xlO~ 5 M0/yr), and hence the accretion luminosity 
is also quite high, L t ~ 5.7 L . The dust temperature increases 
very close to the protostar, but drops down below 20 K within 
a few 1000 AU from the protostar, and then rises towards the 
edge where the core is heated by the ambient radiation field 
(cf. Shirley et al. 2002). Thus, at this early stage, it is neces- 
sary to probe the inner 1000AU to infer the presence of the 
protostar. The peak of the SED has moved to ~ 100 to 120 /mi 
(depending on viewing angle), which implies an effective tem- 
perature T EFF ~ 21 to25K. The dependence of the SED on 
viewing angle is due to the disc around the protostar. 

Time-frame t3, Class object. The bipolar jets from the 
protostar have created hourglass cavities, as can be seen on the 
density plots of Fig. 0] As a consequence, the accretion rate 
onto the protostar decreases to M* ~ 10~ 5 MQ/yr, but its mass 
has increased (to ~ 0.20 Mq), and so the luminosity injected 
into the cloud increases (to L t ~ 27.2 L0). As a result, the 
temperature is > 25 K within 1000 to 1500 AU, and it does 
not fall below 20 K in the inner 5000 AU. The SED peaks at 
~ 80 to 100 /mi (depending on viewing angle), which implies 
r EFF ~ 25 to31K. 

Time-frame t4, Class object. The total luminosity in- 
jected into the cloud is lower (L y ~ 12.3 Lq) due to the re- 
duced accretion rate onto the protostar (M* ~ 2 x 10~ 6 MQ/yr). 
The temperature distribution is similar to the previous time- 
frame (time-frame t2), but slightly cooler due to the reduced 
luminosity. The SED is also similar, but there is a slightly 
stronger dependence on viewing angle. The SED now peaks 
at ~ 90 to 135 /mi, implying T EFF ~ 18 to 28 K. 

Time-frame t5, Class object. Asymmetries in the pattern 
of accretion onto the protostar have given it a small peculiar 
velocity, ~ 0.3 km s _1 , and by t5 it is displaced from the dense 
central region by ~ 100 AU. Consequently the accretion rate 
onto the protostar is reduced (to M+ ~ 4 x lO~ 7 M0/yr) and 
with it the accretion luminosity (to L r ~ 2.5 Lq). As illustrated 
in the temperature cross section parallel to the x — plane on 
Fig.0] the upper half of the core is warmed to > 14 K by the 
displaced protostar. In contrast, the lower half of the core is 
not affected by the protostar, and the temperature here does not 
rise above ~ 10K, apart from in the outer layers of the core, 
which are still heated by the ambient radiation field. As a con- 
sequence of the reduced accretion luminosity, the peak of the 



SED moves back to longer wavelengths (~ 150 /mi) and cooler 
effective temperature (~ 17 K). The displacement of the pro- 
tostar also means that it is now visible at shorter wavelengths 
(1 to 50 /mi), from viewing angles close to the rotation axis. 

5. Summary 

We have simulated the collapse of a turbulent molecular core 
using SPH, and performed 3-dimensional Monte Carlo radia- 
tive transfer simulations at different stages of this collapse, to 
predict dust temperature fields, SEDs and isophotal maps. We 
focus on the initial stages of protostar formation, i.e. just before 
and just after the formation of a protostar, and derive criteria for 
distinguishing between genuine prestellar cores and cores that 
contain very young protostars. As pointed out by Masunaga 
et al. (1998), very young protostars are difficult to observe di- 
rectly, because they are deeply embedded, and may be wrongly 
classified as prestellar cores. Hence it is important to consider 
ways in which their presence might be inferred indirectly. The 
main results of this study are as follows. 

(i) A starless core is heated only by the ambient interstel- 
lar radiation field. Hence its outer regions are warm (~ 20 K) 
and the temperature drops towards the centre of the core (to 
~ 7 K). As the core collapses the central region becomes denser 
and colder (~ 5 K), because the optical depth to the centre in- 
creases. Thus, prestellar cores that are about to form a protostar 
tend to be colder than precollapse prestellar cores. 

(ii) When a protostar first forms in a core, its luminos- 
ity vaporizes the dust in its immediate vicinity, and heats the 
dust that survives up to the dust destruction temperature (~ 
1500to2000K). Initially, only the region within ~ 1000 AU 
of the protostar is hotter than 20 K, but as the collapse proceeds 
the accretion luminosity of the protostar increases and the pres- 
ence of the protostar affects the entire core. At this point, the 
heating of the core is mainly due to the accretion luminosity 
of the protostar, and the effect of the external radiation field is 
secondary. 

(iii) The effective temperature of the SED follows the vari- 
ation of the accretion luminosity, which initially increases due 
to the increase of the protostar mass but then decreases due to 
the decline in the mass accretion rate onto the protostar. Thus 
the SED of a prestellar core peaks at ~ 190 /mi, correspond- 
ing to an effective temperature T EW ~ 13K, but after a pro- 
tostar has formed in the core, the peak of the SED moves to 
shorter wavelengths (~ 80 //m, corresponding to r EFF ~ 31 K) 
as the accretion luminosity increases, and then back to longer 
wavelengths (~ 130 /mi, corresponding to r EFF ~ 19 K) as the 
accretion luminosity decreases. 

(iv) The peak of the SED of a Class object depends only 
weakly on viewing angle. In general the peak is at longer wave- 
length when the Class object is viewed in the plane of the 
surrounding accretion disc. It varies by at most +25 /mi, and is 
always < 150 /mi. Thus, it should be straightforward to distin- 
guish Class objects from prestellar cores, whose SEDs peak 
at ~ 190 /mi. 

(v) In our simulations, a newly-formed protostar cannot 
be observed in the NIR because it is too deeply embedded. 
This result reflects the high densities which our model predicts 
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in the immediate surroundings of a newly-formed protostar. 
Additionally, it depends on the dust opacity there. Thus, the 
presence or absence of NIR emission from very young embed- 
ded protostars (confirmed by the presence of molecular out- 
flows or compact cm radio emission), could be used to con- 
strain the model. 

(vi) At submm and mm wavelengths, isophotal maps of 
cores that contain protostars (Class objects) tend to have 
more structure than maps of prestellar cores, but their emission 
also tends to be more strongly concentrated towards their cen- 
tre. Consequently, after convolving with the beams of current 
telescopes, most of this central structure is lost and cores that 
contain protostars appear almost circular. In contrast, prestellar 
cores tend to have more extended emission and tend to main- 
tain their shape after convolution. 

Based on the above results we propose two criteria for iden- 
tifying cores which - despite appearing to be prestellar - may 
harbour very young protostars: 

(a) They are warm (T > 15 K) as indicated by the peak of 
the SED of the core (/l pea k < 170 /mi). This criterion requires 
that the peak of the SED can be measured to an accuracy of 
~ 30 /mi, which should be feasible with observations in the 
far-IR from ISO (e.g. Ward-Thompson et al. 2002) and the up- 
coming Herschel mission (Andre 2002), and in the mm/submm 
region from SCUBA (Kirk et al. 2005) and IRAM (Motte et al. 
1998). We have presumed that the core is heated by the aver- 
age interstellar radiation field, and the criterion would have to 
be modified for a core irradiated by a stronger or weaker radia- 
tion field. 

(b) Their submm/mm isophotal maps are circular, at least 
in the central 2000 to 4000 AU. 

Using the above criteria, it is relatively easy to identify 
cores that may contain young protostars. These warm, circu- 
lar cores can then be probed by other means (e.g. deep radio 
observations with the VLA, or sensitive NIR observations with 
Spitzer), to ascertain whether they in fact contain very young 
protostars. 
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Fig. 7. 850 jum isophotal maps of 3 different time-frames (first column: 1 1 - collapsing prestellar core; second column: t3 - 
Class object; third column: t5 - Class object), and from three different viewing angles (first row: f? b s = 0°, <p b s = 0°; 
second row: 6> b s = 45°, $> bs = 90°; third row: O b s = 90°, o b s = 0°). Cores with protostars (second and third columns) are 
more centrally condensed than prestellar cores (first column). We note that the axes (x,y) refer to the plane of sky as seen by the 



observer. 
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Fig. 8. As Fig.0 but after convolving with a Gaussian beam having FWHM — 2300 AU, (or equivalently 15" for a core at 150 pc). 
Cores with protostars (middle frame, t3; right-hand frame, t5) appear more circular, but otherwise similar to the prestellar core 
(left-hand frame, tl). 



